Interaction-induced dynamical phase locking of Bose-Einstein condensates 
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We show that interactions result in the emergence of a definite relative-phase between two initially 
incoherent Bose-Einstein condensates. The many-realization interference fringe visibility is universal 
aT - Si2 ~ 1/3 throughout the Josephson interaction regime, as evident from a semiclassical picture. 
Other types of incoherent preparation yield qualitatively different coherence dynamics. 



Ever since their first realization [l[ , interference exper- 
iments between two Bose-Einstein condensates (BECs) 
released from a double-well trap, have raised funda- 
mental questions concerning gauge symmetry breaking 
and the appearance of macroscopic coherence in the 
bose quantum gas. Two extreme cases are usually con- 
trasted: When the condensates are separated, the state 
of the system corresponds to a relative-number squeezed 
state, approaching a Twin-Fock preparation for fully- 
independent and equally populated BECs. The corre- 
sponding interference pattern was predicted to have an 
arbitrary relative phase in each experimental run, vary- 
ing randomly from one realization to another 0, 0, @] • 
Thus, the many-realization average fringe visibility van- 
ishes for such preparation. By contrast, when the two 
condensates are coupled, the initial state corresponds 
to a coherent preparation with a definite relative phase 
between the constituent BECs, reflected by a 'phase- 
locked' reproducible interference pattern and near-unity 
many-realization fringe visibility. Such coherent split- 
ting, maintaining a definite relative phase between the 
condensates, was demonstrated in atom-chip experiments 
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The effect of interactions on the fringe visibility of co- 
herent preparations, has recently attracted much atten- 
tion @, 0, H, & [l(J EH, 0) EH ■ Given time to evolve under 
the influence of inelastic collisions between the atoms, the 
relative-phase of the separated BECs disperses because 
the basis relative-number (Fock) states oscillate with dif- 
ferent frequencies @, 0, S, 0|- Thisprocess has come 
to be known as phase- diffusion [l(| 0j], [l2|, E3- Its dy- 
namics which is closely related to the Josephson effect 
in superconductors, can be studied within the two site 
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Bose-Hubbard model 

H = - JL X + UL 



(1) 



as a function of the characteristic interaction strength 
u = UN/ J, where J and U denote the coupling and 
collisional interaction energies, respectively [1 31 ] . The 
angular momentum operators L x = (a\a2 + a\a\)j2, 



(a[a 2 - a\ax)/{2i), and L z = (hi - h 2 )/2, 



are 



Ly 

defined in terms of bosonic annihilation and creation op- 
erators &i, a\ for particles in condensate % = 1,2, with 
corresponding particle numbers hi = aJSj, satisfying the 



conservation law hi + n2 = N = 21. In the extreme 
strong-interaction Fock regime u > N 2 , single-particle 
coherence is lost on a (U^/l\~ 1 timescale regardless of 
the initial relative phase ip ML Coupling the conden- 
sates results in phase-locking [16j. However, the required 
coupling strength to arrest phase-diffusion depends on ip 
[l3j ]. Whereas a tp — phase is locked already in the 
strong-interaction Josephson regime 1 < u < N 2 , a rela- 
tive phase of ip = tt is only locked in the weak-interaction 
Rabi regime u < 1. 

In this work we consider the effect of interactions on 
the fringe-visibility of the initially incoherent prepara- 
tion. Instead of the initial coherent state obtained by 
fast splitting, it is possible to prepare the relative-number 

We 
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squeezed state by slow separation 
study the buildup of single-particle coherence between 
such separated condensates, due to the combined effect 
of interactions and coupling. In the Fock regime, number 
squeezed states are a good approximation to the system's 
eigenstates so that no coherence may form. However, 
throughout the Josephson regime we find that significant 
coherence may build up, leading to a non-random phase 
distribution in a many-realizations interference experi- 
ment. The resulting fringe visibility attains a uni- 
versal value of ~ 1/3 throughout the Josephson regime, 
in excellent agreement with a semiclassical phase-space 
picture. We also study other, phase-squeezed prepa- 
rations obtained by unitary rotations of the relative- 
number state. Such states are encountered for example, 
in the phase-acquisition stage of Mach-Zendher interfer- 
ometry with number-squeezed inputs [20| . We find that 
fringe visibility buildup for these preparations takes place 
in the Rabi- Josephson transition point and that it is sen- 
sitive to the initial bi-valued relative-phase. 

Below we use for representation the Fock space basis 

states \£, m) a where a — x,y, z, which are the joint eigen- 
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states of L and L ai with I — N/2. We consider the three 



Fock preparations, 

\t,0)a 



N/2 



(2) 



with b x i,2 = (Si ± a 2 )/\/2, b y i,2 = (ai ± ia 2 )/V2, 6 z i,2 = 
ai,2- The states \£,m) x ,y may be obtained by switching 
the coupling and the bias potential between the wells, 



2 



^ 0.2 

"h 

"—'-0.2 
^ 0.6 

<^ 
-0.2 



(a) 



(b) 



(d) 



(e) 



5 10 150 5 10 150 



LUjt 




0.03 



FIG. 1: (Color online) Evolution of fringe visibility with 
iV=1000 particles, starting from the Fock preparations \£, 0) x 
(solid), \£, 0} y (dash-dotted), and \£, 0) z (dashed). The cou- 
pling parameter is u — 0.5 (a), 1 (b) 1.5 (c), 2 (d), 10 (e), and 
100 (f). 



inflicting rapid unitary rotations of the Twin- Fock \£, m) z 
state. 

For all three initial states, the expectation value (L y ), 
vanishes identically throughout the evolution with the 
Hamiltonian ([1]). Thus, the fringe- visibility in a many- 
realizations interference experiment is — \(L x )\/£. 
In Fig. [T] we plot the dynamics of (L x )/£ with N = 1000 
particles, starting from the three Fock preparations, for 
various values of the interaction parameter u. We make 
the following observations: (i) When starting from the 
site number state \£, 0) z interactions lead to the forma- 
tion of a non-vanishing single-particle coherence when 
u < N 2 . This coherence persists well into the Joseph- 
son regime, (ii) Similarly, finite coherence is obtained 
for the \£, 0} x ,y preparations, but only for relatively weak 
interactions u ~ 1. (iii) The coherence evolution for the 
\£, 0) x state exhibits oscillations and beating absent from 
the dynamics of the \£,0) y preparation, approaching a 
different mean value. 

The dynamics of fringe visibility could be interpreted 
using a semiclassical picture 13 , 22|, H^, [U HH . Clasical 
(mean-filed) trajectories are given by the equal energy 
contours of the Gross-Pitaevskii energy functional, 
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sin v cos ip 



(3) 



which is restricted to the unit Bloch sphere. In the Rabi 
regime all motion is essentially linear and the trajecto- 
ries correspond to slightly perturbed Rabi Oscillations. 
In the Josephson regime the spherical phase-space is split 
by a figure-eight separatrix trajectory (Fig[^a)), to a lin- 
ear 'sea' (blue) and two interaction-dominated nonlinear 
'islands' (green). Finally, in the Fock regime the linear 
domain becomes smaller than Planck cell, and therefore 
effectively disappears. In what follows we will assume 
for simplicity that the interaction is repulsive, i.e. u > 0. 
Since u — * — u, (p —f (p + tt implies E — > —E, the phase 
space picture for attractive interactions, is a mirror image 
of the u > contours. 




lJJiI.jI 



^ 

Ly/i ' ' L 3 



o 

E /J 

n 



500 



FIG. 2: (Color online) Schematic illustration of classical 
phase-space energy contours for the three Fock preparations 
l^jO)^ (dashed), \t, 0) y (dash-dotted), and \£,0) z (dotted), are 
shown in (a) with u = 1.5. The corresponding eigenstate 
occupation probabilities P„ are plotted for \£, 0) x (b), 1^,0)^ 
(c), and \£, 0) z (d). Vertical dotted lines in (b)-(d) correspond 
to the separatrix energy. Trajectories of maximum tangency 
(bold) produce caustics in the P n distribution. 



Semiclassical WKB quantization [lj, [23j, |2J, |25| is 
attained by demanding that A(E n ) = (4n/N)(n+ 1/2) 
where A(E) is the phase-space area enclosed by a fixed 
energy E contour, and 4tt/N is the Planck cell. As 
a result, level spacing is determined from the classical 
oscillation frequency w{E) = [A'(E)]^ 1 . The result- 
ing spectrum constitutes in agreement with the classi- 
cal phase-space structure: (a) low-energy sea levels ex- 
tending from E_ = —£J with Josephson-frequency spac- 
ing ojj = w(EJ) = y/(J + NU)J, (b) separatrix lev- 
els around E x — £J with spacing w x = 2ujj/ log(7V 2 /u), 
and (c) high-energy island levels, approaching E+ w £ 2 U, 
with characteristic spacing lo+ w NU between nearly de- 
generate pairs. We note parenthetically that in the strict 
classical limit {N — > 00 keeping u fixed) the separatrix 
frequency vanishes [2(| l27l| bu t due to the slow log(TV) 
convergence to classicality [ll| , for finite N it only differs 
from the Josephson frequency by a logarithmic factor. In 
the Rabi regime, we have E + < E x and the entire spec- 
trum consists of sea levels, whereas in the Fock regime 
the ground state's energy lies above E x so that all levels 
reside in the islands. 

To understand the dynamics of the three Fock prepa- 
rations, it is useful to consider their Wigner function 
p(8,ip). The Wigner functions depicting the \£,0) a 
preparations, lie along the corresponding great circles 
around the a axes, as shown in Fig. EUa). The Wigner 
functions of the BHH eigenstates are concentrated 
along the contour lines of the Gross-Pitaevskii classical 
energy functional. Thus, the expansion coefficients of 
the \£, Q) a preparation, in terms of the BHH eigenstates 
\E n ) can be estimated by the semiclassical prescription 
P{E n ) = \(E n \iP)\ 2 = traced pW). 

In Fig. [2fb)-(d) we plot the coefficients P{E n ) for the 
\£,fyx,y,z preparations, for u = 1.5 in the Josephson 
regime, obtained by direct numerical diagonalization of 
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FIG. 4: (Color online) Acquired time-averaged coherence as a 
function of the interaction parameter u for the Fock prepara- 
tions \i, Q) x (solid), \l, 0) y (dash-dotted), and \£, 0) z (dashed). 



FIG. 3: (Color online) Frequency distributions P(uj) for the 
preparations \(.,0) a , a — x (top), y (middle), and z (bottom) 
with iV=1000 particles, and u = 0.5 (a), 1.5 (b), and 10 (c). 
Dotted lines mark the Josephson frequency uij whereas the 
low frequency limit corresponds to the separatrix frequency 
cj x , logarithmically approaching zero with increasing N. 



the BHH. It is clear from the phase-space landscape of 
Fig. [21a) that the eigenstate expansion should be quite 
different for the three initial states. The state \£, 0} x lies 
entirely in the sea, overlapping a narrow band of lin- 
ear eigenstates. By contrast, the state \£, 0) y straddles 
the entire spectrum, including sea, separatrix, and is- 
land levels, whereas the \£, 0) z state consists of all sea 
levels up to the separatrix energy, but does not project 
at all onto the nonlinear islands. This observation is re- 
flected well in the numerical results with the expected 
extent of the local density of states and caustics ob- 
tained for trajectories of tangency with the initial Wigncr 
distribution. For example, for the \£, 0) x preparation, 
the highest expansion coefficients in Fig. [2jb) are ob- 
tained for the energy contours E n » E(ir/2,ir/2) = 
and E n » E(n,ip) = uNJ/A = 375 J (for u = 1.5 and 
N = 1000), marked by bold (blue) lines in Fig.EJa). 

In order to relate the probability distributions of 
Fig. [2jb)-(d) to the time evolution depicted in Fig. [TJ 
we plot the frequency distribution P{uj n ) with oj n = 
E n+ i — E n , for the three initial Fock states, in Fig [3j In 
this picture, frequencies extend from the separatrix lev- 
els with to ~ v x , through the equally-spaced to ~ toj low 
energy sea levels, to the maximal level spacing uj ~ lo + 
at the top of the nonlinear islands. In the Rabi regime 
(FigE^a)), the \£, 0) x preparation corresponds to a nar- 
row distribution of sea levels with two dominant frequen- 
cies corresponding to the classical trajectories tangen- 
tial to the L x — great circle. Hence beating is ob- 
served around the Josephson frequency. In comparison, 
the states 1^,0)^ have at their disposal the entire sea 



frequency range, allowing for the continuous buildup of 
single particle coherence. 

In the Josephson regime (Fig G2Jb),(c)) the separatrix 
and islands appear. The states \£, 0) x and \£, 0) y begin to 
penetrate the islands and project into the high-frequency 
regime at u — 2 and u = 1, respectively. Consequently, 
if u is large enough, (L x ) never attains a significant mag- 
nitude (solid and dash-dotted lines in Fig. |4j. By con- 
trast, the frequency-span of the preparation \£, 0} z re- 
mains within the w x < lu < luj range and time-averaged 
coherence of ~ 0.35 — 0.40 at a relative phase ip = 
is maintained throughout the Josephson regime (Fig. [T] 
and dashed line in Fig. [4|. Finally, in the Fock regime, 
the entire phase space is nonlinear, the \£, 0) z preparation 
becomes a narrow distribution of (island) eigenstates of 
the BHH, and the interaction does not induce dynamical 
phase-locking. 

The value of the interaction-induced fringe-visibility 
obtained for initially separated condensates can be de- 
duced from the classical time evolution of their Wigner 
distribution p(9, ip). The Wigner distribution of the twin 
Fock state \£,0) z (Fig. E^a)) lies along the L z — great 
circle, passing throughout the linear sea region of phase 
space, up to the separatrix trajectory. Due to the vari- 
ation of characteristic frequency from tu x near ip = tt to 
ujj near p = 0, propagation in time results in a 'spi- 
ral' motion of the Wigner distribution (Fig. [Spa)), which 
for ujt ^> 1 becomes spread out throughout the sea 
(Fig. E^c)). The phase-distribution, obtained by tracing 
p{9, p) over 9 is not uniform since all classical sea tra- 
jectories pass through ip = whereas only the separatrix 
trajectory passes through ip = ir. The average value of 
the coherence (L x ) for this nearly-uniform distribution is 
given by 



J* cos [9 x (p)} cos(p)dtp ^ 
J* cos [0x0)] (<p)dip 



(4) 
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FIG. 5: (Color online) Popagation of 1001 classical trajecto- 
ries with u = 100, starting from initial conditions correspond- 
ing to the Wigner distribution of the twin-Fock preparaion: 
(a) t = 0, (b) uijt = 2, (c) ujt = 50. Gray lines correspond 
to the mean-field energy contours of Eq. @ . The mean value 
of L x over all points (solid) is plotted in (d) and compared to 
the numerical quantum dynamics of Fig. [ljf ) (dashed). 

where x (ip) is the separatrix energy contour, 

- cos 2 [9 x (ip)} - sin [0 X (<^)] cos ^ = 1 . 

In the Joscphson regime we can approximate 
cos[9 x ((p)] = + cos ip)/u. Substituting this 

approximate separatrix line into Eq. ((4]) , we obtain that 
(L x )ujjt^>i — This value is in good agreement with 
the numerical results of Fig. [4] Averaging numerically 
over the classical distribution, we find that the coher- 
ence dynamics overlaps the numerically-exact quantum 
calculation (Fig. E^d)), with a quasi-equilibrium value of 
(L x )u>jt^>i = 0-37£ for u — 100. This value is universal 
throughout the Josephson regime, regardless of the exact 
values of u and N. 

In conclusion, starting from a relative-number state 
of equally populated BECs (a Twin-Fock state), inter- 
actions in the Josephson regime 1 < u < N 2 , result in 
the build-up of single-particle coherence. This phase- 
locking process is robust, insensitive to the exact value 
of the coupling, and has the opposite effect compared 
to phase-diffusion. Significant average fringe visibility 
°f 9i2 ~ 1/3 is generically obtained throughout the 
Josephson regime, regardless of interaction parameters 
or particle number , in excellent agreement with semi- 
classical analysis. Interactions thus serve to select a 
non-random relative-phase in the weak merging of ini- 
tially fully-separated condensates. The proposed mech- 
anism is fundamentally different from the deterministic, 
single-condensate reflection fringes obtained in collisions 
of effectively immiscible BECs [28[ and requires a much 



weaker interaction strength, typical of current BEC in- 
terference experiments. For the phase-squeezed states 
obtained by unitary rotations of the Twin-Fock state, we 
find phase-sensitive dynamics of the fringe-visibility in 
the weak interaction regime, also explained to excellent 
accuracy within the semiclassical picture. 
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